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ABSTRACT 

The observed UV continuum slope of star forming galaxies is strongly affected by the 
presence of dust. Its observation is then a potentially valuable diagnostic of dust atten- 
uation, particularly at high-redshift where other diagnostics are currently inaccesible. 
Interpreting the observed UV continuum slope in the context of dust attenuation is of- 
ten achieved assuming the empirically calibrated Meurer et al. (1999) relation. Implicit 
in this relation is the assumption of an intrinsic UV continuum slope (j3 — —2.23). 
However, results from numerical simulations suggest that the intrinsic UV continuum 
slopes of high-redshift star forming galaxies are bluer than this, and moreover vary 
with redshift. Using values of the intrinsic slope predicted by numerical models of 
galaxy formation combined with a Calzetti et al. (2000) reddening law we infer UV 
attenuations (Ai 500 ) 0.35 — 0.5 mag (Ay: 0.14 — 0.2 mag assuming Calzetti et al. 2000) 
greater than simply assuming the Meurer relation. This has significant implications 
for the inferred amount of dust attenuation at very-high (z ps 7) redshift given cur- 
rent observational constraints on j3, combined with the Meurer relation, suggest dust 
attenuation to be virtually zero in all but the most luminous systems. 

Key words: galaxies: evolution - galaxies: formation - galaxies: starburst - galaxies: 
high-redshift - ultraviolet: galaxies 
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INTRODUCTION 



The intrinsic UV luminosity of a galaxy is a frequently used 
diagnostic of star formation (e.g. Wilkins et al. 2012b). How- 
ever, UV emission is easily attenuated by dust. Determin- 
ing the level of dust attenuation is then critical to accu- 
rately derive intrinsic UV luminosities and thus star forma- 
tion rates. At high-redshift, far-IR observations, which probe 
dust-reprocessed UV emission (and thus obscured star for- 
mation), are currently inaccesible for all but the most lumi- 
nous galaxies. 

Using observations taken with Wide Field Camera 3 
(WFC3) on the Hubble Space Telescope means it is now rou- 
tinely possible to identify star forming galaxies at z — 7 — 8 
(e.g. Oesch et al. 2010, Bouwens et al. 2010b, Bunker et al. 
2010, Wilkins et al. 2010, Finkelstein et al. 2010, McLure et 
al. 2010, Wilkins et al. 2011a, Lorenzoni et al. 2011, Bouwens 
et al. 2011, McLure et al. 2011, Finkelstein et al. 2012b, 
Lorenzoni et al. 2012) and even to z ~ 10 (e.g. Oesch et 



al. 2012). These deep near-IR observations also allow us to 
probe the rest-frame UV continuum (with at least 2 filters 
unattenuated by the Lyman-cv break) at z = 4 — 7 allowing 
the measurement of the UV continuum slope f3 (definecrl 
such that f\ oc A' 3 ). 

The observed UV continuum slope has been demon- 
strated to provide a reasonable diagnostic of moderate dust 
attenuation in local star forming galaxies (e.g. Meurer et 
al. 1999, hereafter M99). By virtue of being UV selected, 
the UV continuum slopes of high-redshift galaxies are rela- 
tively easy to measure, even for the faintest sources, as they 
require observations in only a single additional filter rela- 
tive to the three required to robustly identify high-redshift 
galaxies. Consequently the UV continuum slope is one of the 
few accessible diagnostics of galaxies that are available at 
very high-redshift and has been subject to significant inter- 
est (e.g. Stanway et al. 2005, Bouwens et al. 2009, Bunker et 
al. 2010, Bouwens et al. 2010a, Wilkins et al. 2011b, Dunlop 
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et al. 2012a, Bouwens et al. 2012, Finkelstein et al. 2012a, 
Rogers, McLure & Dunlop 2012, Dunlop et al. 2012b). 

The observed UV continuum slope (of star forming 
galaxies) can be interpreted, in the context of dust, with 
knowledge of both the intrinsic UV continuum slope (of star 
forming galaxies) and the reddening law. A shortcut is to use 
a relation empirically derived using observations of the UV 
continuum slope and far-IR emission such as that proposed 
by M99 (though this implicitly assumes an intrinsic UV con- 
tinuum slope). However, the M99 relation is calibrated at 
low-redshift, where the properties of star forming galaxies, 
and in particular the intrinsic UV continuum slopes may be 
significantly different from the systems we can now study at 
high-redshift. Indeed recent predictions from the GALFORM 
semi-analytical model of galaxy formation (see Wilkins et al. 
2012a) suggest the intrinsic slope varies with redshift and is 
bluer than that implicit in the Meurer relation. 

In this work we employ a large cosmological hydrody- 
namic simulation (MassiveBlack-Il) to make predictions for 
the intrinsic UV continuum slope in high-redshift star form- 
ing galaxies and contrast the attenuation inferred assuming 
these with those from the M99 relation. This work is or- 
ganised as follows: In Section [2] we describe how the rela- 
tionship between the inferred level of dust attenuation and 
the observed UV continuum slope is affected by various as- 



sumptions including the choice of intrinsic slope ({2.1 1 and 
reddening law ({2.2 1. In Section[3]we consider the implica- 
tions of changing the these assumptions on the level of dust 
attenuation inferred from recent observations of the UV con- 
tinuum slopes of high-redshift galaxies. Finally in Section [1] 
we present our conclusions. 

Throughout this work magnitudes are calculated using 
the AB system (Oke & Gunn 1983). We assume Salpeter 
(1955) stellar initial mass function (IMF), i.e.: £(m) = 
diV/dm oc m" 2 - 35 . 



2 RELATING THE OBSERVED UV 
CONTINUUM SLOPE TO DUST 
ATTENUATION 

The effect of applying either an SMC-Zi&e or Calzetti et al. 
(2000) law to the UV continuum causes the observed slope 
/3 bs to redden relative to the intrinsic value. Using the ob- 
served UV continuum slope as a diagnostic of dust attenu- 
ation in high-redshift star forming galaxies then essentially 
requires knowledge of two factors: (1) the intrinsic UV con- 
tinuum slope, and (2) the reddening law. 

The relationship between the observed slope and the 
attenuation at a given wavelength A can then be defined 
simply as: 



A x = D x x [/3 obs - /3 int ], 



(1) 



where /3j nt is the intrinsic UV continuum slope, and /3 b s is 
the observed slope. The factor Dx relates a change in /3 to 
the attenuation at A (Ax), i.e. Dx = dAx/d/3, and is sensitive 
to the choice of reddening law. In the following sections we 
discuss each of these factors in turn. 



2.1 Factors affecting the intrinsic UV continuum 
slope and predictions from numerical models 
of galaxy formation 

The intrinsic UV continuum slope of a stellar population is 
determined by the joint distribution of stellar masses, ages 
and metallicities. These properties are, in turn, determined 
by the initial mass function (IMF) together with the star for- 
mation and metal enrichment histories. The effect of each of 
these factors was explored in detail by Wilkins et al. (2012a). 



2.1.1 Predictions from numerical models of galaxy 
formation 

By using a galaxy formation model to predict both the 
star formation and metal enrichment histories we can pre- 
dict the intrinsic UV continuum slope of galaxies (see also 
Wilkins et al. 2012a). Here we make use of a state-of-the- 
art cosmological hydrodynamic simulation of structure for- 
mation: MassiveBlack-11. The MassiveBlack-11 simulation is 
performed using the cosmological TreePM-Smooth Particle 
Hydrodynamics (SPH) code P-Gadget, a hybrid version 
of the parallel code Gadget2 (Springel 2005) tailored to 
run on the new generation of Petaflop scale supercomput- 
ers. MassiveBlack-11 includes N pa , r = 2 x 1792 3 « 11.5 bil- 
lion particles in a volume of 10 6 Mpc 3 //i 3 (100Mpc//i on a 
side), is run to z = 0, and includes not only gravity and 
hydrodynamics but also additional physics for star forma- 
tion (Springel & Hernquist 2003), black holes and associated 
feedback processes (Di Matteo et al. 2008, Di Matteo et al. 
2012). Spectral energy distributions (SEDs) of galaxies are 
generated using the stellar population synthesis (SPS) code 
Pegase.2 (Fioc & Rocca-Volmerange 1997, 1999). This sim- 
ulation has been demonstrated to reproduce the observed 
intrinsic UV luminosity function (LF) at z = 5 — 8 (with the 
exception of the faint-end slope at z = 5, see Wilkins et al. 
in-prep) . 

In Figure [T] the median intrinsic UV continuum slope, 
as a function of the intrinsic UV luminosity, for galaxies z £ 
{5, 6, 7, 8, 9, 10} predicted by this simulation are shown. For 
simplicity the UV continuum slope /3 is measured using two 
artificial rest-frame UV bandpasses uniformly encompassing 
1300 — 2700A. This range is similar to the range accesible 
for observations of galaxies at z ~ 3 — 6 using WFC3 on 
the Hubble Space Telescope. However, at present, for sources 
located at z ~ 7 and above a smaller range of the rest-frame 
UV continuum is accesible using WFC3 observations. With 
this in mind we also determine the UV continuum slope of 
galaxies at z ~ 7 using a shorter rest-frame UV wavelength 
baseline 1300 - 2 100 A roughly corresponding to the range 
probed by the WFC3 J/i25to and .ff/ieoio filters at z = 7. 
The effect of changing to this reduced wavelength baseline 
reduces the intrinsic value of /3 by ~ 0.03. For clarity Figure 
[TJpresents the UV continuum slope over the full 1300- 2700 A 
baseline, however when we come to interpret observations at 
z ~ 7 we make use of the predicted intrinsic values of j3 using 
the smaller baseline. This reduced wavelength baseline can 
also have a significant effect when determining the level of 
dust attenuation inferred from an observed slope through a 
change to the parameter Dx - This is discussed in more detail 
in p^l 

Apparent in Figure[T]is a strong trend with redshift such 
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Figure 1. The intrinsic UV continuum slope /3- lnt as a function 
of UV luminosity predicted by the MassiveBlack-Il simulation. 
The upper panel shows both a 2D number density histogram and 
the median slope (and 16 th — 84 th percentile range) in several 
luminosity bins. The shade of the 2D density histogram denotes 
the number of galaxies contributing to each bin on a linear scale. 
Where there are fewer than 25 galaxies contributing to each bin 
on the 2D histogram the galaxies are plotted individually. The 
lower panel shows only the median slope as a function of absolute 
magnitude but for a range of redshifts (z £ {5, 6, 7, 8, 9, 10}) high- 
lighting the important redshift evolution. The horizontal dotted 
line in both panels denotes the intrinsic slope implicit in the M99 
relation (/3 = -2.23). 



that galaxies at higher-redshift typically have bluer UV con- 
tinuum slopes (Pint « —2.55 at z = 8 while /3i nt ~ —2.40 
at z = 5). This reflects the typically younger stellar pop- 
ulations (and thus the presence of more high-mass stars) 
and lower metallicities in higher-redshift systems (Wilkins 
et al. in prep). There is also a trend with intrinsic UV lumi- 
nosity such that more luminous objects typically have red- 
der intrinsic UV continuum slopes. This is again a conse- 
quence of older stellar populations and higher metallicities, 
though this time in more luminous systems. However, this 
effect is fairly small with /3- m t increasing by less than 0.1 
as M1500 = — 18 — > — 21 at z — 7. These predictions closely 
match the luminosity and redshift trends predicted using the 
GALFORM semi-analytical galaxy formation model utilised in 
Wilkins et al. (2012a). However, the normalisation differs 
slightly due to a difference in choice of initial mass function 
(Wilkins et al. 2012a assumed a Kennicutt 1983 IMF). 



2.2 The choice of reddening law 

The factor, D\, relating the change in /3 due to the UV at- 
tenuation can be easily calculated for reddening laws in the 
literature. For the Calzetti et al. (2000) attenuation curve 
we find this factor to be fi600 = 2.10. For the Large and 
Small Magellanic Clouds (LMC and SMC) extinction curves 
(using Pei et al. 1992) we obtain values of Diqoo = 1.63 
and -D1600 = 1.12 respectively while for the core collapse 
supernovae (CCSNe) curve of Bianchi & Schneider (2007) 
we find D1600 = 2.49. The implication of the smaller val- 
ues of D\ obtained for the SMC and LMC laws is that for 
the same observed UV continuum slope a smaller (and sig- 
nificantly smaller assuming the LMC law) UV attenuation 
would be inferred. These values are based on the same wave- 
length baseline used to measure the UV continuum slope 
(1300 — 2700A). If the luminosity baseline is reduced to 
(1300-2100A), as is appropriate when interpreting the UV 
continuum slope of sources observed by WFC3 at z ~ 7, 
these values become 1.84, 0.96, 1.47 and 1.90 for the Calzetti 
et al. (2000), SMC, LMC, and CCSNe curves respectively. 
The strong sensitivity of Dieoo f° r the CCSNe curve reflects 
the structure in the curve around 2000 - 2500A. 



2.3 The Meurer relation 

An alternative to explicitly assuming an intrinsic slope and 
reddening law is to use an empirical relation based on obser- 
vations of the UV continuum slope and some observational 
measurement of the UV attenuation (typically from far-IR 
observations). A prominent and widely used example of such 
an empirical calibration is the Meurer et al. (1999) relation. 
This relation is defined as: 



Ai 



4.43 + 1.99 x /3 t 



(2) 



Re-arranging this relation to match the form of Eqn. [T] re- 
veals the M99 relation implicitly assumes an intrinsic slope 
of /3 bs = —2.23 and factor relating the change in /3 to 
UV attenuation of D1600 = 1.99. The value of -D1600 is un- 
surprisingly similar to that derived assuming the Calzetti 
et al. (2000) reddening law. The implicit intrinsic slope 
(/Sobs = —2.23) however is significantly redder than that 
predicted for star forming galaxies at high-redshift. Consid- 
ering the similar values of D1600 this then suggests assuming 
the intrinsic slope(s) predicted by MassiveBlack-11 combined 
with the Calzetti et al. (2000) reddening law would yield 
larger UV attenuations for the same observed value of p. 



3 IMPLICATIONS FOR INFERRING DUST 
ATTENUATION AT HIGH-REDSHIFT 

In the previous section we described how the observed UV 
continuum slope is sensitive to the choice of intrinsic slope 
and reddening law. We now investigate how these different 
assumptions affect the interpretation of observations of the 
UV continuum slope at high-redshift. 



3.1 The observed UV continuum slope 

As noted in the introduction numerous groups have recently 
studied the observed UV continuum slopes of high-redshift 
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star forming (Lyman-break selectecQ galaxies (e.g. Stanway 
et al. 2005, Bouwens et al. 2009, Bunker et al. 2010, Bouwens 
et al. 2010a, Finkelstein et al. 2010, Wilkins et al. 2011a, 
Dunlop et al. 2012a, Finkelstein et al. 2012a, Bouwens et al. 
2012, Rogers et al. 2012, Dunlop et al. 2012b). A compilation 
of measurements of f3 from some of these studies (at z = 
5 — 7) is shown in Figure [2] 

The measurements from different groups are largely 
consistent, with typical observed slopes of /3 ~ — 2. There are 
however continuing conflicting claims regarding a correla- 
tion with redshift and/or luminosity. For example, Bouwens 
et al. (2009), Wilkins et al. (2011b), Bouwens et al. (2012), 
and Finkelstein et al. (2012) all find evidence of a correla- 
tion with redshift (such that galaxies at lower redshift have 
redder observed slopes) while Dunlop et al. (2012a) do not. 
Bouwens et al. (2009), Bouwens et al. (2010), Wilkins et al. 
(2011b), and Bouwens et al. (2012) find a significant cor- 
relation with luminosity while Dunlop et al. (2012a) and 
Dunlop et al. (2012b) do not. Finkelstein et al. (2012a) do 
not observe a strong relation between the observed /3 and 
UV luminosity at z ~ 4 or 5 in their full sampl^] though 
they do note a trend at z ~ 7 (though this may be the result 
of selection bias). Several possible explanations, relating to 
both the methodology used to measure /3, the UV luminos- 
itjj^] and the bias in the sample are suggested to reconcile 
this apparent discrepancy in both Bouwens et al. (2012) and 
Finkelstein et al. (2012a). 
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3.2 The effect of changing assumptions on the 
inferred level of dust attenuation 

While there is some disagreement between observational 
studies in this paper we are simply interested in exploring 
the effect of different assumptions on the inferred level of 
dust attenuation. With this in mind we consider only the 
Bouwens et al. (2012) observations in our subsequent analy- 
sis, though note it would be equally valid to use the Finkel- 
stein et al. (2012a) results considering the similar sample 
size and luminosity baseline. 

We then consider four sets of assumptions controlling 
the conversion of the observed UV continuum slope to a UV 
attenuation: (1) the M99 relation, (2) the intrinsic slope(s) 
predicted by MassiveBlack-II combined with a Calzetti et al. 
(2000) reddening law (C00+MBII scenario), (3) the intrinsic 
slope(s) predicted by MassiveBlack-II combined with a Pei 
et al. (1992) empirical SMC reddening law (SMC+MBII), 
and (4) the intrinsic slope(s) predicted by MassiveBlack- 
II combined with collapse supernovae (CCSNe) curve of 
Bianchi & Schneider (2007) (CCSNe+MBII). 

The level of (UV) attenuation (j4isoo) inferred from the 
Bouwens et al. (2012) observations using these four sets of 

2 While several studies use a photometric redshift selection tech- 
nique most of the significance comes from the strong Lyman-a 
break feature thus they are de facto LBGs. 

3 However, Finkelstein et al. (2012a) note if they restrict the 
faintest luminosity bin to include only galaxies from the deeper 
HUDF dataset they do find evidence for a correlation with lumi- 
nosity. 

4 Finkelstein et al. (2012a) do find a correlation with stellar mass, 
this is found both observationally and predicted by MassiveBlack- 
II to correlate with UV luminosity. 



Figure 2. Measurements of the observed UV continuum slope in 
the literature at z £ {5, 6, 7}. For the measurements from Dunlop 
et al. (2012b) we utilise their determination of /3 using a power- 
law methodology (i.e. column 4 of table 1 of Dunlop et al. 2012b). 
Finkelstein et al. (2012a) present their result in three luminosity 
bins referenced on the local value of L* (L < 0.25xL*, 0.25xL* < 
L < 0.75 X L* , L > 0.75 x L*). The median is represented by the 
solid line while the shaded regions show the uncertainties. For the 
L < 0.25 x L* and L > 0.75 x L* bins the minimum and maximum 
bounds respectively are chosen to reflect the approximate range 
of contributing galaxy luminosities. 



assumptions are shown in both Figure [3] (where we show the 
inferred attenuation as a function of observed UV luminosity 
for the three redshifts bins) and Figure [4] (where we show 
the inferred attenuation as a function redshift for several 
luminosity bins). As noted previously for the z ~ 7 sample 
we utilise the intrinsic UV continuum and values of -D1500 
measured assuming a smaller wavelength baseline given the 
coverage of the WFC3 observations. 

The difference in the inferred attenuation (A1500) be- 
tween the M99 relation and C00+MB11 scenario is largely 
fixed as a function of luminosity reflecting the similar func- 
tional form of the Calzetti et al. (2000) reddening curve and 
that implicit in the M99 relation (i.e. the value of D\). This 
difference does however evolve with redshift due to evolution 
of the intrinsic slope predicted from the simulations (and the 
non-evolution of the intrinsic slope implicit in the M99 re- 
lation) . At 2 ~ 5 the difference in A1500 is ^ 0.35 mag and 
increases to ~ 0.5 mag at z ~ 7. Given that the applica- 
tion of the M99 relation to observations at z fa 7 results in 
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zero (or even negative) attenuations at fainter luminosities 
this difference has a strong impact. This increase in the in- 
ferred attenuation brings the level of attenuation even for 
the faintest luminosities consistent with zero. The effect of 
assuming a CCSNe curve is also similar to the COO+MBII 
scenario considering the similar values of Disoo- 

The difference between assuming the M99 relation and 
the SMC+MBII scenario is more nuanced. While the value 
of -D1500 is larger for the Calzetti law than the SMC law 
(i.e. for a similar increase in f3 a greater attenuation is in- 
ferred using the Calzetti law) the intrinsic slope predicted 
by MBII is bluer than that implicit in the M99 relation. For 
very-blue observed values of /3 (< —2) this difference in the 
implicit intrinsic slope is the most important factor and re- 
sults in a greater level of attenuation inferred assuming the 
SMC+MBII scenario than utilising M99. A consequence of 
this is that at low luminosities the level of dust inferred 
assuming the SMC+MBII scenario does not drastically fall 
below zero. At redder observed values of /3 (> —2) this trend 
reverses and the difference in D1500 dominates. Thus at very- 
high luminosities the level of attenuation inferred assuming 
the SMC+MBII scenario is lower than simply assuming M99 
(though this is only true for the most luminous point at 
z = 5-6). 

As can be seen in Figure |4] the adoption of a decreasing 
intrinsic UV continuum slope to higher redshift (as opposed 
to the fixed value implicit in the M99 relation) also weakens, 
though does not destroy, the inferred (decreasing) trend of 
UV attenuation with redshift. 
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4 CONCLUSIONS 

It has been established, at low-redshift, that the UV con- 
tinuum slope is a useful diagnostic of dust attenuation in 
typical star forming galaxies. However, the relation between 
the level of attenuation and the observed slope is sensitive 
to both the assumed intrinsic slope and reddening law. 

While it is possible to use an empirically determined 
calibration (such as the Meurer et al. 1999 relation) this 
nevertheless implicitly assumes an intrinsic UV continuum 
slope (in this case /3i nt = —2.23) which may not be valid for 
very-high redshift where the UV continuum slope is the only 
accesible diagnostic of dust attenuation. Indeed, we find that 
results from the large cosmological hydrodynamical simula- 
tion MassiveBlack-Il suggest the intrinsic slope both varies 
strongly with redshift, and, at z > 4 is bluer than the value 
implicit in the Meurer relation. 

Making use of the intrinsic UV continuum slopes pre- 
dicted by MassiveBlack-II (and the Calzetti et al. 2000 red- 
dening law) we re-interpret recent observations of the UV 
continuum slope finding levels of dust attenuation 0.35 — 
0.5 mag (where the difference is sensitive to the redshift and 
luminosity) higher than simply applying the Meurer rela- 
tion. If instead we utilise an SMC-like redding law (com- 
bined with the intrinsic slope predicted by MassiveBlack-W) 
we find levels of attenuation similar to those inferred using 
the Meurer relation with the exception of the lowest lumi- 
nosities (where we infer a higher level of attenuation). 



Figure 3. UV attenuation inferred from observations (Bouwcns 
et al. 2012: B12) of the UV continuum slope at z e {5,6,7} 
using three different assumptions: (1) The empirically determined 
M99 relation, (2) intrinsic UV continuum slopes predicted from 
MassiveBlack-II and a Calzetti et al. (2000) reddening law, (3) 
intrinsic UV continuum slopes predicted from MassiveBlack-II 
and a SMC type reddening law, and (4) intrinsic UV continuum 
slopes predicted from MassiveBlack-II and the CCSNe reddening 
law of Bianchi & Schneider (2007). 
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